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We show that Abe's general pseudoadditivity for entropy prescribed by 

thermal equilibrium in nonextensive systems holds not only for entropy, but 

also for energy. The application of this general pseudoadditivity to Tsallis 

entropy tells us that the factorization of the probability of a composite system 

into product of the probabilities of the subsystems is just a consequence of 

the existence of thermal equilibrium and not due to the independence of the 

subsystems. 
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Suppose an isolated system composed of two subsystems 1 and 2. One of the basic 
assumptions of thermodynamics is the existence of an equilibrium state in the compound 
system at which T\ = T 2 , where T is the absolute temperature. This so called zeroth 
law is obeyed by Boltzmann statistical mechanics having additive entropy and energy. For 
nonextensive statistical mechanics (NSM) |T[, the validity of this law is a more subtle affair 
which depends on the relationships aSu = fi2(aSi,aS 2 ) and bEi 2 = gi 2 (bEi,bE 2 ), 
where S is entropy and E the internal energy. The numbers in the indexes of the functions 
indicate the dependence of the latter, a and b are constants used to make each variable 
dimensionless in the above equalities. From now on, we let a = b = 1 so that S and E 
become dimensionless. 

The zeroth law has been established, in a approximate way 0-§J, for NSM with Tsallis 
entropy and additive energy E 12 = Ex + E 2 . Recently, Abe JFj furthermore find a general 
pseudoadditivity of entropy required by the existence of thermal equilibrium for additive 
energy : 

H(S 12 ) = H(S X ) + H{S 2 ) + \ S H(S 1 )H(S 2 ), (1) 

where H is certain different iable function satisfying H(0) =0 and As a constant depending 
on the nature of S. For Tsallis entropy H can be proved to be the identity function. 
Eq.(|l]) is very interesting because it can be considered as a general criterion of pertinent 
nonextensive entropies for equilibrium systems and may help to understand Tsallis nonex- 
tensive statistical mechanics and to find other nonextensive thermostatistics obeying the 
zeroth law. 

Nonextensive entropy is a consequence of long range correlations or complex (fractal or 
chaotic) space-time. When interactions are no more limited between or on the walls of the 
containers of subsystems, energy may be nonadditive. This issue has been widely discussed 
and energy nonadditivity was clearly shown for some cases with long range interactions 
|. So, inspired by Abe's work, we naturally ask the following question : what is the kind 



of energy nonadditivity that satisfies the requirement of existence of thermal equilibrium? In 



this letter, along the line of Abe, we show that the above equilibrium pseudoadditivity holds 
not only for entropy, but also for internal energy in more general cases where energy-type 
thermodynamic variables are not extensive. 

To establish the zeroth law, we make a small variation of the total entropy dS\2 given by 

dS 12 = ^^dS x + ^rdS 2} (2) 

and a variation of total energy E12 given by 

At thermal equilibrium, dS\2 = should hold and dEyi = results from energy conservation, 
leading to 

9/12 dSl -dE x = - d -i^^dE 2 , (4) 



and 



dSi dE 1 8S 2 dE 2 



dE; dEl = 'dE 2 dE2 - (5) 



So we obtain 



d/12 ac d/i2 

(6) 
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Equilibrium means that this above equation yields the zeroth law which is in general given 
by following equality : 

Ft = F 2 (7) 

where F x and F 2 are the same functions depending on subsystem- 1 and —2, respectively 
This constraint from the zeroth law implies following factorizations of the derivatives in 
Eq.® : 

= 012^1^2, (8) 



-Qg- = 012^1^2, (9) 
<%12 



dEx 

and 



Ol2fJ>lV2, (10) 



9 X2 vxH2- (11) 



dE 2 

where 0, u, u, 9, v and fi are certain functions of the subsystems indicated by the indexes. 
From Eq.(|) we can write 

dSi dS 2 

^dE~-^dE 2 (12) 
with £i = an d £2 — • This shows that Eq.(||) to (|TT|) are really the most general 
forms of the derivatives of / and g satisfying Eq.(^). 

We find that all the calculations of Abe [7j hold for / as well as for g. So we can replace 
S by E in Eq. ([!]). For a system containing iV subsystems in equilibrium, we have 

JV 

ln[l + X x H x (x 12 .. N )} = $>[1 + \ x H x ( Xi )}. (13) 

i=i 

where x can be entropy S or energy E. In general, Hs and As are different from He and 
A^, respectively. 

We mention here as example two nonextensive cases with Tsallis entropy (S = — 1 '^ r q pq , 
p is density operator) where the zeroth law is claimed to be verified. 

1. In the Tsallis nonextensive statistics with escort probability, Hs and He are the 
identity function, \ s = (I — q) and Xe = |l|fJ-Pl (i.e. E is extensive and 
S 12 = S 1 + S 2 + (l-q)S 1 S 2 ). 

2. Another possible case is with Tsallis entropy combined with a so called incomplete 
normalization p~7^[T9[] where H x is the identity function and X x = q — 1 for both 
entropy and energy (i.e. x\ 2 = x\ + x 2 + (q — l)a;ia;2 or, according to Eq.(|13|), xi 2 ..n = 
[rL=i( 1 + A3:X ') jj with N subsystems. In this case, the zeroth law can hold without 
approximation, making it possible to establish an exact nonextensive thermodynamics. 



An example of systems satisfying pseudoadditivity of entropy is given by Abe [|7| with 
H x (x) = ^fx for black hole entropy proportional to its horizon area. This discussion should 



hold for energy as well according to the first law of thermodynamics for black hole |20| if 
electromagnetic work is absent. So it would indeed be interesting to study black hole within 
a generalized thermodynamics with nonadditive entropy and energy as well in view of the 



difficulty with Boltzmann statistics due to the presence of thermal (infrared) divergence [21 



Another example can be given with the long-range ferromagnetic spin model of which the 
internal energy is given by E(N,T) = c{T)N N x "' d l i^pl where N is the number of 



-a/d 

particles in the model supposed additive, d the dimension of space, a the exponent in the 
factor l/r a of the long range potential and c(T) certain function of temperature T. 

When N — » oo and a > d, E(N,T) = c(T)N is additive. On the other hand, when a = d, 
E(N,T) = c(T)N\nN. If we put, e.g., H E (E) = {exp[e E / c ^ N ] - e}/e and = 1, then 
E satisfies Eq.(|) or Eq.flT3|). When < a < d, E(N, T) = c{T)N 2 - a / d , we can choose 
H E {E) = [^} 1/{2 ~ a/d) and \ E = for energy to satisfy Eq.(|) or Eq.(|T3|). There exist 
other choices, e.g. H E (E) = ew{[^\ l/{2 ' a/d) } ~ 1 and \ E = 1. 

As a matter of fact, it seems that, for a given explicit relation between energy or en- 
tropy and number of subsystems or volume supposed additive (i.e. V\ 2 = V\ + V 2 and 
N12 = N\ + N2), the finding of a function H satisfying Eq.(|l|) is a trivial affair. The most 
essential contribution of Abe's work, in our opinion, is that Eq.(|l]) finally makes it clear 
that the factorization of the compound probability of a composite system into product of 
the probabilities of the subsystems is a consequence of the existence of thermodynamic equi- 
librium if Tsallis entropy applies, because H is identity function here. It is straightforward 
to show that : S 12 = S 1 + S 2 + {1- q)S 1 S 2 = -±=£& - + (1 - q) 1 -^^ 1 -^ 



1-<J l-q 1 V" ^) 1-q 1-q 

i-Tr( PiP2 )« _ i-Trp^ 2 



( _ - i-q 12 ' which means p\ 2 = p\p 2 [or that p?-(l,2) = p^(l)p^(2) implies 
Pij(l, 2) = pi(l)pj(2) where Pi(l) or Pj(2) is the probability for subsystem-1 or -2 to be at 
state i or j and Py(l, 2) is the probability for the composite system to be at the product state 
So the product probability is rather due to thermal equilibrium instead of independence of 
noninteracting or weakly interacting subsystems as claimed or insinuated in most of the cur- 



rent publications. On this basis, we can finally free ourselves from the paradox of addressing 
noninteracting independent systems with nonadditive entropy due to long range interactions. 
Additive energy based on independent systems with only short range interactions does not 
conform with the spirit of nonextensive statistical mechanics. 

The law of product probability must be respected for any thermodynamic system in equi- 
librium. This means that, for NSM with the Tsallis g-exponential distribution, nonextensive 
energy is needed for exact treatments of nonextensive interacting systems. In this sense, all 
treatments within NSM with additive energy should be viewed as a kind of extensive approx- 
imation and should be proceeded with great care because they give sometimes very different 
results from the treatment respecting probability factorization or product state fl9| , |2~2"|| . This 
above "equilibrium interpretation" of the factorization hypothesis of compound probability 
may have important consequences on the applications of NSM to many-body systems. More 



detailed discussions on this issue are given in others papers of ours [23 



Summing up, we have applied Abe's method of finding general entropy pseudoadditivity 
to a more general case where both entropy and energy are nonextensive. Under thermal 
equilibrium, the energy of a nonextensive composite system also obeys Abe-type pseudoad- 
ditivity. 

We would like to thank Professor Sumiyoshi Abe for critical reading of this manuscript 
and important comments. Thanks are also due to Professor S. Ruffo for valuable discussions. 
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